Statistics of excitons in quantum dots and the resulting microcavity emission spectra 
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I. INTRODUCTION 

Semiconductor quantum dots (QDs)i are a leading 
technology for the investigation of the quantum realm. 
They offer exciting possibilities for quantum computa- 
tion and are important candidates for the next genera- 
tion of light emitters. In most cases, the best control 
of the states of the confined carriers in QDs is obtained 
through coupling to light 2 . This light- matter interaction 
can be considerably enhanced by including the dot in 
a microcavity, with pillars^, photonic crystals^ and mi- 
crodisks^ bein g the currently favoured realizations. Ref- 
erences f^ l ll r J| describe the first reports, in each of these 
structures, of vacuum field Rabi splitting, whereby one 
excitation is transferred back and forth between the light 
and the matter fields. This contrasts with the weak cou- 
pling regime previously studied^*i where only quantita- 
tive perturbations of the dynamics occur, such as reduc- 
tions in the lifetimes of the dot excitations (Purcell ef- 
fect). In the case of strong coupling, however, the coher- 
ent exchange of energy merges the light and matter ex- 
citations into a new entity. This is commonly referred to 
as an exciton-polariton in semiconductor physics^ with 
an important example being the two-dimensional polari- 
tons in planar microcavities, first observed by Weisbuch 
et al£ In cavity quantum electrodynamics (cQED), the 
equivalent concept is the dressed state of atoms by the 
quantised electromagnetic field. 

In QDs, optical interband excitations create electron- 
hole pairs or excitons, confined by a three-dimensional 
potential which makes their energy spectrum discrete. 
If this potential is much stronger than the bulk exciton 
binding energy, and if the size of the dot is smaller than 
the corresponding exciton Bohr radius, the Coulomb in- 
teraction between electrons and holes can be considered 
as a perturbation. For the lowest exciton states, this is 
the fermionic limit where the Pauli exclusion principle 
dominates. In the opposite limit, if the confining poten- 
tial is weak or the size of the dot is much greater than the 
exciton Bohr radius, the exciton is quantized as a whole 
particle. In this case, the bosonic nature of excitons is 



expected to prevail over fermionic nature of individual 
electrons and holes. An important question for the de- 
scription of emission from QDs embedded into cavities 
in the strong coupling regime is whether the dot excita- 
tions coupled to light behave like fcrmions or like bosons. 
Here we address the question of which statistics (Bosc- 
Einstein, Fermi-Dirac or a variation thereof) best de- 
scribes excitons in QDs. This is a question which is very 
topical in view of the recent experimental achievements, 
and which has elicited substantial theoretical works in 
the past, in connection with the possibility of exciton 
Bose condensation. 

In this paper we derive the exciton creation operator 
in a QD which allows the calculation of nonlinear optical 
spectra of QDs in microcavities. The model we develop 
takes into account the saturation of the transition due to 
Pauli exclusion alone and does not attempt to solve the 
complex manybody problem which arises when Coulomb 
interactions between excitons are included. Hence the 
model is most accurate in describing the departure from 
ideal bosonic behaviour in large dots rather than near 
the fermionic limit in small dots. We analyze the dot 
size effect on the statistics of excitons and demonstrate 
the transition from the fermionic to bosonic regime. To 
motivate this, we begin by summarizing how the cou- 
pling of light modes with fermionic and bosonic material 
excitations differ. 

The Rabi doublet, with splitting amplitude 2Kg as 
shown on Fig. QJa), is well accounted for theoretically 
by the coupling with strength fig of two quantized oscil- 
lators a and b both obeying Bose algebra, 



[a, a 1 "] = 1, 



(1) 



and equivalently for b. We shall describe this well known 
and elementary case in detail as it provides the founda- 
tion for most of what follows. Neglecting off-resonant 
terms like aW, the hamiltonian readsiS 



H = Tiuj(a)a + tfb) + fig(ab i + a f 6) 



(2) 



We assumed degenerate energies fiuj for the two oscilla- 
tors, which will not affect our qualitative results, while 
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simplifying considerably the analytical expressions. One 
oscillator, say a, describes the light field while the other, 
b, describes a bosonic matter field. The analysis of (0 
can be made directly in the bare state basis \i,j) with i 
excitations in the matter field and j in the photonic field, 
i, j G N. The value of this approach is that the ex- 
citation, loss and dephasing processes generally pertain 
to the bare particles. For instance matter excitations 
are usually created by an external source (pumping) and 
light excitations can be lost by transmission through the 
cavity mirror. This physics is best expressed in the bare 
states basis. 

On the other hand (0 assumes a straightforward ex- 
pression in the basis of so-called dressed states which di- 
agonalises the Hamiltonian to read: 



H = (Hlu — hg)p^p + (Hu> + hg)q'q 



(3) 



where p^ and create a coherent superposition of bare 
states, respectively in and out of phase: 



p = (a + b)/V2, 



(a-b)/V2. 



(4) 



For clarity we shall note the dressed states, i.e., 

the eigenstates of Q with i dressed particles of en- 
ergy Tiuj — hg and j of energy hu> + hg. The manifold 
of states with n excitations thus reads, in the dressed 
states basis: 
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Its energy diagram appears on the left of Fig. ^ for man- 
ifolds with zero (vacuum), one, two and seven excita- 
tions. When an excitation escapes the system while in 
manifold 7i„, a transition is made to the neighbouring 
manifold TL n -i and the energy difference is carried away, 
cither by the leaking out of a cavity photon, or through 
exciton emission into a radiative mode other than that of 
the cavity, or a non-radiative process. The detailed anal- 
ysis of such processes requires a dynamical study, but 
as the cavity mode radiation spectra can be computed 
with the knowledge of only the energy level diagrams, we 
shall keep our analysis to this level for the present work. 
The important feature of this dissipation is that, though 
such processes involve a or & (rather than p or q), they 
nevertheless still result in removing one excitation out of 
one of the oscillators. Hence only transitions from 
to \\i — 1, j) or — 1) are allowed, bringing away, re- 
spectively, fiui + hg and huj — hg of energy, accounting for 
the so-called Rabi doublet (provided the initial i and j 
are nonzero in which case only one transition is allowed) . 
From the algebraic point of view, this of course follows 
straightforwardly from Q and orthogonality of the basis 
states. Physically it comes from the fact that, as in the 
classical case, the coupled system acts as two indepen- 
dent oscillators vibrating with frequencies u> ± g. In the 
case of vacuum field Rabi splitting, a single excitation is 
shared between the two fields, and so the manifold Hi 
is connected to the single line of the vacuum manifold. 
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FIG. 1: (Color online) Energy diagrams of the two limiting 
cases of dressed bosons (left) and fermions (right). In the first 
case the nth manifold has constant energy splitting of 2hg 
between all states and couples to the (n — l)th manifold by 
removal of a quantum of excitation with energy hu±hg which 
leads to the Rabi doublet, Fig. (a), with splitting 2hg. In the 
second case, each manifold is two-fold with a splitting which 
increases like a square root. All four transitions are allowed, 
leading to the Mollow triplet, Fig. (b), for high values of n 
when the two middle transitions are close in energy. The 
distance from the central peak goes like 2hg^/n and the ratio 
of peaks is 1 : 2. The two lowest manifolds (in blue) are 
the same in both cases, making vacuum field Rabi splitting 
insensitive to the statistics. 
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In this case there is obviously no possibility beyond a 
doublet. 

Different physics occurs when the excitations of the 
material are described by fermionic rather than bosonic 
statistics. In the case of cavity QED the material is usu- 
ally a beam of atoms passing through the cavity, and 
a single excitation is the independent responses of the 
atoms to the light field excitation. The simplest sit- 
uation is that of a dilute atomic beam where a single 
atom (driven at resonance so that it appears as a two- 
levels system) is coupled to a Fock state of light with a 
large number of photons. This case is described theoret- 
ically by the Jaynes-Cummings modelii, in which a (the 
radiation field) remains a Bose operator but b becomes 
a fermionic operator which describes two-levels systems, 

— > <j + with : 

ff + = (o o) ■ (6) 

Then the atom must be in either the ground or the 
excited state, allowing for the manifolds 

W n = {|0,n),|l,n-1)} (7) 

provided that n > 1 . The associated energy diagrams ap- 
pear on the right of Fig. with two states in each man- 
ifold (in our conventions |0, n) refers to the bare states 
with atom in ground state and n photons, while |1, n — 1) 
has atom in excited state and n—1 photons). For the res- 
onant condition, where the two-levels transition matches 
the cavity photon energy, the dressed states for this man- 
ifold are split by an energy ^fnhg. In the general case, 
all four transitions between the states in manifolds TL n 
and Ti n -i arc possible, which results in a quadruplet. It 
is hard to resolve this quadruplet, but it has been done 
in Fourier transform of time resolved experiments^. It 
is simpler to consider photoluminescence directly under 
continuous excitation at high intensity (where the fluctu- 
ations of particles number have little effect). In this case, 
with the two intermediate energies are almost de- 

generate and a triplet is obtained with its central peak 
being about twice as high as the two satellites. This is 
the Mollow triplet of resonance fluorescenceii 

Thus we arc faced with two limiting cases, one is a pure 
bosonic limit with equally spaced dressed states resulting 
in the linear Rabi doublet, the other the pure fermionic 
limit with pairs of dressed states of increasing splitting 
within a manifold but decreasing energy difference be- 
tween two successive manifolds, giving rise to the Mol- 
low triplet. In many of the strong coupling experiments 
conducted so far, and in all the reports concerning semi- 
conductor QDs, only one single excitation is exchanged 
coherently, so that the states are dressed by the vacuum 
of the electromagnetic field, resulting in the Rabi dou- 
blet. However, at this level of excitation, there is com- 
plete agreement between the bosonic and fermionic mod- 
els, with both providing a good description of the experi- 
mental observations. The prospect of stronger pumping, 



with more than one excitation shared between the two 
fields, makes it important to understand whether a re- 
alistic semiconductor QD will correspond to the bosonic 
or fermionic case, or something intermediate between the 
two. 

We review and discuss some of the more significant 
achievements in this field in section [n] In section IIIII 
we lay down a general formalism for building the exci- 
ton creation operator. In section Hvl we study two limit- 
ing cases which resemble Bose and Fermi statistics. We 
show how in the general case the luminescence behaviour 
interpolates between these two limits which we have al- 
ready discussed, and calculate the second order coher- 
ence of the emitted light. In section [V] we draw the 
experimental consequences of the various statistics and 
discuss how the spectra obtained allow a qualitative un- 
derstanding of how excitations distribute themselves in 
the excitonic field. In the final section we conclude and 
discuss briefly ongoing work to refine the modelling of 
the excitonic wavefunction. 



II. EXCITONS AS QUASI-PARTICLES 

The generic optical excitation in an intrinsic semicon- 
ductor is the electron-hole pair. In bulk, the two op- 
positely charged particles can be strongly correlated by 
the Coulomb interaction and bound as hydrogenic states 
(Wannier-Mott exciton). Although finding binding en- 
ergies and wavefunctions for the single exciton case is 
a difficult problem in various geometries as far 

as the vacuum coupling limit is concerned, the exciton 
field operator which links the two manifolds Ti\ and vac- 
uum always assumes the simple form a+ regardless of the 
details of the exciton. As a particle constituted of two 
fermions, the exciton is commonly regarded as a boson, 
from consideration of the angular momentum addition 
rules and the spin-statistics theorem. For a single par- 
ticle, this is an exact statement, albeit a trivial one. It 
is however generally agreed to hold at small densities^. 
In this case the Rabi doublet is obtained, as observed 
experimentallyiiiiii. 

At higher excitation power, the problem assumes con- 
siderable complexity as well as fundamental importance 
for physical applications. Already at the next higher 
excitation — with one more electron, hole or electron-hole 
pair added to the first exciton — the situation offers rich 
and various phenomena both in weaki&iS and strongSS 
coupling, owing to the underlying complexity of the ex- 
citonic states. In this work we shall be concerned with 
resonant optical pumping, so that excitations are created 
in pairs and the system always remains electrically neu- 
tral. We shall describe as an exciton any state of an 
electron-hole pair, whether it is an atom-like Is state or 
has both particles independently quantised^ and in a 
more general sense we shall also use the term for any 
combination of particles which takes part in the exci- 
tonic phase. Indeed the excitonic phase with more than 
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one pair requires at the most accurate level a description 
in terms of an excitonic complex, e.g., in term of biexci- 
tons/bipolaritons for strong coupling of two electrons and 
two holes with lights. This comes from the Coulomb in- 
teraction which links all charge carriers together and, in 
a most fundamental way, also from the antisymmetry of 
the wavcfunction which demands a sign change whenever 
two identical fcrmions (electrons or holes) are swapped in 
the system. However in some configurations, especially in 
planar cavities, a widely accepted hypothesis of bosonic 
behavior of excitons and the derived polaritons has been 
investigated for effects such as the exciton boser— , po- 
lariton amplifier 22 i 23 i 24 and polariton lasero 25 ! 26 . The in- 
ternal structure of the exciton which gives rise to both 
deviations from the Bose-statistics and interactions of the 
electron-hole pairs is then expressed as an effective repul- 
sive force in a bosonised Hamiltonian (due to Coulomb 
interaction and Pauli effect in the form of phase-space 
filling or exchange interaction^LS). 

This bosonic approach for excitons met early oppo- 
sition in favour of an analysis in the electron-hole ba- 
s j s 29 i 30 4 3^ Combcscot and co-workers investigated the 
possibility of bosonisation of excitons^SiSi^iS and con- 
cluded against it. They point out its internal incon- 
sistency, as the same interaction binds the underlying 
fermions, and therefore defines the exciton, while also 
being responsible for exciton-exciton scattering; this is 
inconsistent with the indistinguishability of the particles. 
These authors introduce the "proteon" as the paradigm 
for Bose-like composite particle and propose a formal- 
ism ("commutation technics") which essentially relies 
on evaluating quantum correlators in the fermion basis 
with operators linked through the single exciton wave- 
function. The importance of Fermi statistics of the un- 
derlying constituents has also been pointed out by Rom- 
bouts and co-workers in connection to atomic and exci- 
tonic condensates^^. In both cases the composite "bo- 
son" creation operator reads 

^=E^l k (8) 

k 

in term of Uk, Qc the fermion operators 

{ck, er£,} = 5 k ,k', {ok, ov} = (9) 

(same for st), respectively for an electron and hole (or 
proton in the atomic case^l) of momentum k. This op- 
erator creates excitons with the center-of-mass momen- 
tum K = in a system with translational invariance. It 
was appreciated long ago^£ that the operators and B 
are no-longer exact bosonic operators, but instead satisfy, 

[B,B*] = 1-D, (10) 

where D is a nonzero operator, though with small matrix 
elements at low exciton densities. 

Our approach is based on a definition similar to Eq. (JSJ 
for the exciton operator. Instead of analysing the devia- 
tions in the commutation relationship i|l(J|) we derive the 



matrix elements of the operator B^ . The direct anal- 
ysis of these matrix elements allows us to trace depar- 
tures from Bose-statistics and to investigate the transi- 
tions between bosonic and fermionic behaviours of exci- 
tons i n QDs. We shall c ompare our results with those of 
Refs. |32l33l34l35l36l37j . For convenience, we will refer to 
these two sets of publications through the names of their 
first authors, keeping in mind that they are co-authored 
papers, as listed in the references. 

III. FORMALISM 

We consider a QD, which localises the excitation in 
real space. Thus our main departure from Combescot 
and Rombouts is that our exciton creation operator, X', 
is expressed in real space and without the zero center-of- 
mass momentum restriction, 

n e ,nu 

where cr„ c and <; nh are fermion creation operators, cf. 10, 
for an electron and a hole in state \<p% ) and ¥>£ h ), re- 
spectively: 

<|o) = K), 4jo> = |<>, (12) 

with |0) denoting both the electron and hole vacuum 
fields. We carry out the analysis in real space with the 
set of basis wavefunctions 

^(r e ) = (r e KJ and <(r h ) = falO (13) 

with r e and the positions of the electron and hole, 
respectively. Subscripts n e and rih are multi-indices enu- 
merating all quantum numbers of electrons and holes. 
The specifics of the three-dimensional confinement man- 
ifests itself in the discrete character of n e and n% compo- 
nents. 

We restrict our considerations to the direct band semi- 
conductor with non-degenerate valence band. Such a sit- 
uation can be experimentally achieved in QDs formed 
in conventional III-V or II- VI semiconductors, where the 
light-hole levels lie far, in energy, heavy-hole ones due 
to the effects of strain and size quantization along the 
growth axis 2 !. Therefore, only electron-heavy hole exci- 
tons need to be considered. Moreover, we will neglect 
the spin degree of freedom of the electron-hole pair and 
assume all carriers to be spin polarized. This can be 
realized by pumping the system with light of definite cir- 
cular polarization, as spin-lattice relaxation is known to 
be very inefficient in QDs. 

The (single) exciton wavefunction \tp) results from the 
application of X' on the vacuum. In real space coordi- 
nates: 

(r e ,r h \tp) = (p{r e ,r h ) = ^ C n ^ nh ip e ne {r e )^ h {v h ). 

n e ,nh 

(14) 
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At this stage we do not specify the wavefunction (that 
is, the coefficients C„ ej „ h ), which depends on various fac- 
tors such as the dot geometry, electron and hole effective 
masses and dielectric constant. Rather, we consider the n 
excitons state which results from successive excitation of 
the system through X^: 



I*, 



(Xt)»|0>, 



(15) 



We later discuss in more details what approximations are 
being made here. For now we proceed by normalising this 
wavefunction 



where, by definition of the normalization constant 



Mn = V(*n|*n>. 



(16) 



(17) 



The creation operator can now be obtained explic- 
itly. We call a n the non-zero matrix element which lies 
below the diagonal in the excitons representation: 

a n = (n\X*\n-l), (18) 

which, by comparing Eqs. (|15H18|) turns out to be 



(19) 



We now undertake to link a n with the coeffi- 
cients C netnh , which assume a specific form only when 
the system itself has been characterised. The general re- 
lationship is more easily obtained in the real space than 
with the operator representation (|llfl . Indeed the non- 
normalized n-excitons wavefunction assumes the simple 
form of a Slater determinant 



*n(r ei , . . .,T en ,r hl ,. ■ .,Th n ) 



det [^(re^rfc)], 

l<ij<n 



(20) 

explicitly ensuring the antisymmetry of \ff n upon ex- 
change of two identical fermions (holes and electrons), 
as results from the anticommutation rule ® in the nth 
power of operator X\ cf. I|llfl . 

The determinant can be computed explicitly, by expan- 
sion of its minors which results in the recurrent relation 



-. n m— 1 



m—l 



3=0 



with Ao = 1 and (3 m the irreducible m-excitons overlap 
integrals, 1 < m < n: 



m—l 

IK 

i=l 



(T ei ,r hi )<p(r ei ,r hi+1 ) (p*(r £m ,r hm )p(r em ,r hl ) dr ei . ..dr em dr hl . ..dr hri 



(22) 



The determinant can also be solved by direct combinato- 
rial evaluation, counting all combinations which can be 
factored out as products of (3 m . The result expressed in 
this way reads: 

p(n) n 

^ = E^)IIfi W ( 23 ) 

r\— 1 m—l 

where p(n) is the partition function of n (number of ways 
to write n as a sum of positive integers, i.e., as an integer 
partition of n) and v-qii) is the number of times that i £ 
N appears in the 77 th partition of n. The coefficients C n 
read: 

cm '-'nf- 3 ? 1 **)^ 

where the product is taken over the integers i%i which 
enter in the 77 th partition of n, i.e., n = X)i n i- 

The procedure to calculate the matrix elements of 
the creation operator is as follows: One starts from the 



envelope function y>(r e ,r^) for a single exciton. Then 
one calculates all overlap integrals j3 m as given by 12211 . 
for 1 < to < n where n is the highest manifold to be 
accessed. Then the norms can be computed, the more 
practical way being recursively with (|21J) . Finally the 
matrix elements a„ are obtained as the successive norms 
ratio, cf. i|18|) . Once a n are known the emission spectra 
can be calculated with ease. We note here that the nu- 
merical computation of the j3 m and a n values needs to 
be carried out with great care. The cancellation of the 
large numbers of terms involved in Eq. (|21|l requires a 
very high-precision evaluation of /3 m . 

Although derived from the formulation in real 
space H20fl . the recurrent relation l|21|l . or its analytical 
solution (|23|1 . is a property of fermion pairs, so it applies 
to the fermionic operators in (|SJ) as well. The core of the 
mathematical results contained in these two expressions 
has in fact been obtained by Combesco^ through di- 
rect evaluation with the operator algebra involved in B' . 
The only difference with her approach and ours is that 
her corresponding quantity p m (which she notes a m ) ap- 
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pears as a series in the (reciprocal space) wavefunction 



k 



(24) 



as opposed to the overlap integral Q22|). 

With such a simple expression as (|24p. Combescot ei aZ. 
have been able to obtain approximate analytical forms for 
Is states of excitons in both 3D and 2D. The sum over 
reciprocal space is approximated as an integral in the 
continuum limit, and since 4>k depends on k = |k| only, 
this becomes £ k |0k| 2m -> / V k \(j)k\ 2m dk with V k = ink 2 
or 2irk, respectively, yielding^ 



^-t 16 (8m-5)! 



(8m- 2)!! 



2y r 



3m — 1 



in 3D, (25) 
in 2D, (26) 



with z and y dimensionless parameters involving the ratio 
of Bohr radius ob to the size of the system. In Ref. [3^ , 
this continuous approximation is used even when the sys- 
tem size becomes so small that the quantization becomes 
noticeable, as a result of which a negative norm is ob- 
tained even for a well defined wavefunction. However, 
the authors interpret incorrectly this negative norm as 
the result of an unphysical exciton wavefunction and thus 
the breakdown of the bosonic picture. To demonstrate 
this, we have computed numerically the values of a m by 
direct summation of Eq. (|24fl and used them to evaluate 
normalization constants Af n . Some results of the compu- 
tation are shown in Fig. |2 The inset shows the depen- 
dence of A/i = o"i vs. the ratio of crystal size L to 3D 
exciton Bohr radius ae- The continuous approximation 
Eq. (|25|l gives <j\ = 1, shown by the dashed line in the 
inset. The solid line displays the value of <j\ computed 
exactly with (|24|l : it approaches 1 when L/ob — > oo but 
departs strongly from 1 as the size of the system reduces, 
since the real wavefunction changes as the system shrinks 
about it. At the very least this needs to be taken into ac- 
count by correcting the normalization constant of ^>k ; as 
otherwise unphysical result may arise. The main part of 
the figure shows the two-excitons normalization A/2 cal- 
culated by both methods. We find that for all L/clb, the 
exact A/2 never becomes negative in clear constrast with 
the behavior obtained with the approximated value 1)25( 1 . 
Hence it is not possible to use the appearance of a nega- 
tive norm as a criterion for bosonic breakdown. 



IV. LIMITING CASES 

We shall not attempt in this paper to go through the 
lenghty and complicated task of the numerical calculation 
of the exciton creation operator matrix elements for a re- 
alistic QD. Rather, we consider here a model wavefunc- 
tion which can be integrated analytically and illustrates 




FIG. 2: Normalization constant A/2 vs. ratio of crystal size 
L to 3D exciton Bohr radius hb- Dashed line shows the re- 
sult using the approximated expression i12ol for <j m and solid 
line using the exact evaluation of the sum on the mesh in 
reciprocal space. The negative value of the norm used as a 
criterion for boson wavefunction breakdown in Ref. [32I is an 
artifact of this approximation. With the genuine wavefunc- 
tion, our procedure yields a non-arbitrary way to consider how 
the wavefunction vanishing norm affects the bosonic charac- 
ter of the excitation. Inset displays A/i vs. Lja-g, computed 
on mesh (solid) and according to Eq. 125H (dashes). 



some expected typical behaviours. Before this, we com- 
ment briefly on the limiting cases of Bose-Einstein and 
Fermi-Dirac statistics and how they can be recovered in 
the general setting of this section. 

This will be made most clear through consideration 
of the explicit case of two excitons (n = 2). Then the 
wavefunction reads 

* 2 (r ei ,r ea ,r fcl ,r Aa ) = ip(r ei ,r hl )(p(r e2 ,r h2 ) 

- p(r ei ,r h2 )tp(r e2 ,r hl ) (27) 

with its normalization constant 1(17)1 readily obtained as 

■W? = J l*2(r ei ,r e2 ,r/ ll ,r, l2 )| 2 dr ei . ..dr h2 = 2 - 2f3 2 

(28) 

where /3 2j the two-excitons overlap integral, reads explic- 
itly 



02 = J <p(r ei ,T hl )<p(T es ,r h2 ) 

ip(r ei ,r, l2 )tp(r e2 ,r hl )dr ei ...dr h2 . (29) 

This integral is the signature of the composite nature of 
the exciton. The minus sign in 1)28)) results from the Pauli 
principle: two fermions (electrons and holes) cannot oc- 
cupy the same state. Assuming tp(r e ,Th) is normalised, 
A/i = 1, so according to l(19)l. 



a 2 



V2 - 2/3 2 



(30) 
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Since < /?2 < 1 this is smaller than or equal to s/2, the 
corresponding matrix element of a true boson creation 
operator. This result has a transparent physical meaning: 
since two identical fermions from two excitons cannot be 
in the same quantum state, it is "harder" to create two 
real excitons, where underlying structure is probed, than 
two ideal bosons. We note that if L is the QD lateral 
dimension, /3a ~ {as/h) 2 <C 1 when L ^> ob- Thus in 
large QDs the overlap of excitonic wavefunctions is small, 
so ct2 ~ y/2 and the bosonic limit is recovered. On the 
other hand, in a small QD, where Coulomb interaction 
is unimportant compared to the dot potential confining 
the carriers, the electron and hole can be regarded as 
quantized separately: 



<p(r e ,r h ) = ip e {Y e )y h (r h ) 



(31) 



In this case all j3 m = 1 and subsequently all a m = 
at the exception of a\ = 1. This is the fermionic limit 
where maps to the Pauli matrix er + , cf. (jfJJ). 

We now turn to the general case of arbitrary sized 
QDs, interpolating between the (small) fermionic and 
(large) bosonic limits. We assume a Gaussian form for 
the wavefunction which allows to evaluate analytically 
all the required quantities. As numerical accuracy is not 
the chief goal of this work we further assume in-plane 
coordinates x and y to be uncorrelated to ease the com- 
putations. The wavefunction reads: 



tp(r e ,r h ) = Ccxp(-7 e r 2 
properly normalized with 



C 



leh 



(32) 



(33) 



provided that j eh G [—2y/%rfa,0] with ~/ e , lh > 0. The 
7 parameters allow to interpolate between the large and 
small dot limits within the same wavefunction. To con- 
nect these parameters j e , 7/1 and 7^ to physical quan- 
tities, (|32|l is regarded as a trial wavefunction which is 
to minimize the hamiltonian Hqd confining the electron 
and hole in a quadratic potential where they interact 
through Coulomb interaction^: 



H, 



QD 



^ \ 2m; 2 



? 9 



(34) 



Here Pi the momentum operator for the electron and 
hole, i = e, h, respectively, m e , rrih the electron and hole 
masses, uj the frequency which characterises the strength 
of the confining potential, e the charge of the electron 
and e the background dielectric constant screening the 
Coulomb interaction. This hamiltonian defines the two 
length scales of our problem, the 2D Bohr radius 0,3 and 
the dot size L: 



,h 2 



«B 



L 




(35a) 



(35b) 




FIG. 3: Ratio of parameters — 7^ and 7 (with 7 = 7 e = 
7h) as a function of L/cib- For large dots where L > ab, 
— 7eh/7 ~ 2 which corresponds to the bosonic limit where 
the electron and hole are strongly correlated. For shallow 
dots where L <g ae, — 7eh/7 ~ with electron and hole 
quantized separately. The transition is shown as the result 
of a variational procedure, with an abrupt transition when 
the dot size becomes comparable to the Bohr Radius. 



where /i = m e mh / \m e + mu) is the reduced mass of the 
electron-hole pair. To simplify the following discussion 
we assume that m e = rrih, resulting in "f e = 7;, = 7. 
The trial wavefunction l|32() separates as ip(r e ,rh) = 
C$(R)</>(r) where r = r e — Yh is the radius-vector of rel- 
ative motion and R = (r e + r^)/2 is the center-of-mass 
position: 



$ (R) = ^2(27 + le h ) exp ( _ R2[27 + 7efe]) (36a) 



, \ V 2 7 - leh I 2 

(r) = j== — exp I — r 



/27T 



27 - leh 



(36b) 



Eq. (|36a|l is an eigenstate of the center-of-mass energy 
operator and equating its parameters with those of the 
exact solution yields the relationship 27 + j e h = 2/L 2 . 
This constrain allows to minimise (|36b|) with respect to 
a single parameter, a = —~f e h/2 + 1/(2L 2 ), which even- 
tually amounts to minimise 4ae/a 2 + a^a 2 / L 4 — /a. 
Doing so we have obtained the ratio —^ e h/l as a function 
of L /ob displayed on Fig. [3] The transition from bosonic 
to fermionic regime is seen to occur sharply when the dot 
size becomes commensurable with the Bohr radius. For 
large dots, i.e., for large values of L/clb, the ratio is well 
approximated by the expression 



- 7eh /7 = 2 - {a B /Lf 



(37) 



so that in the limit a-^/L —* 0, Eq. I|32|) reads ip(r e ,Th) oc 
cxp(— (y^%r e — ^/T/Tr/j) 2 ) with vanishing normalization 
constant. This mimics a free exciton in an infinite quan- 
tum well. It corresponds to the bosonic case. On the 
other hand, if L is small compared to the Bohr ra- 
dius, with 7 e h — > 0, the limit l|31|l is recovered with 
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tp oc exp(— 7 e Tg) cxp(— 7^r?). This corresponds to the 
fermionic case. 

One can readily check that l|32[l gives, in the case 
7 e /j — > — 2^/7 e 7/j, an cxciton binding energy which is 
smaller by only 20% than that calculated with a hydro- 
genic wavcfunction, which shows that the Gaussian ap- 
proximation should be tolerable for qualitative and semi- 
quantitative results. Moreover, its form corresponds to 
the general shape of a trial wavcfunction for an exciton 
in an arbitrary QDjS with the only difference that we 
take the Gaussian expression instead of a Bohr exponen- 
tial for the wavefunction of the relative motion of the 
electron and hole. This compromise to numerical accu- 
racy allows on the other hand to obtain analytical ex- 
pressions for all the key parameters, starting with the 
overlap integrals (|22|l which take a simple form in terms 
of multivariate Gaussians: 



and 



ftn 



where 



C 2m J exp(-x T Ax)dx J cxp(-y T Ay) dy (38) 



X (*^ei ) *^e2 ) • * * ; *^e m j ^hi j ^h,2 i • • • j ^h m ) (39tl) 

y T = (ve 1 ,y e2 , ■ ■ ■ ,Ve m ,Vhi,yh 2 , ■ ■ -,yh m ) (39b) 

are the 2m dimensional vectors which encapsulate all the 
degrees of freedom of the m excitons-complex, and A is 
a positive definite symmetric matrix which equates 11' 21) 
and l|38fl . i.e., which satisfies 



2 m 



X T Ax = 27 e E x i + 2>lh E X i + "fehX m X m+ i 
i—1 i=m-\-l 

rn rri—1 

+ leh E X i X m+i + leh E X i X ™+i+l (40) 



i=l 



i=l 



and likewise for y (to simplify notation we have not writ- 
ten an index m on x, y and A, but these naturally scale 
with /3 m ). The identity for 2m- fold Gaussian integrals 



cxp(— x T Ax) <ix 



VdetA 



(41) 



allows us to obtain an analytical expression for /3 m , al- 
though a cumbersome one. The determinant of the ma- 
trix A reads 



m m — k / \ i 

dctA = W n V V(-i) Lm/2J+ U m (fc,o (2^ ) 



k=0 1=0 

Here we introduced a quantity 



(42) 



Am(k,l) = A' m (k,l) 

j n 

+ E (^n-riK I - i) - I ~ 0) (43) 



A' m (k,l) 



p(0 

E 

rj=l 



m — I 



(44) 



with fc g]0, m], ^ G [0,m] and p(Z) and already in- 

troduced as the partition function of I and the number 
of occurence of i in its rjth partition. For the case k — 
the finite size of the matrix implies a special rule which 
reads A m (0, 1) = 4<5 mi ;<5 m = 2 ,o- Together with Ipffil. i|4*T|) 
and l|42|l . expression l|38[l provides the /3 m in the Gaussian 
approximation. One can see the considerable complex- 
ity of the expressions despite the simplicity of the model 
wavcfunction. Once again, even a numerical treatment 
meets with difficulties owing to manipulations of series 
of large quantities which sum to small values. We had to 
turn to exact algebraic computations to obtain a n coef- 
ficients free from numerical artifacts. Before we present 
the numerical results, we once again turn to the lim- 
iting case of a large QD (L 3> ob) where the bosonic 
behaviour may be expected, putting again for simplic- 
ity 7 e = 7^ = 7. When Eq. I|37[) holds, rather lengthy 
algebraic manipulations yield 



ft, 



(2a B /L) 2(m - 1) 



in- 



(45) 



This approximation is valid for small values of ob /L and 
for m -C L/a-Q. In computations of a n with n -C £/ob, 
the denominator in Eq. I|45|l plays a minor role and 



vWl-2(n-l)(a B /£) 5 



(46) 



Thus in the small-77, limit the matrix elements of the 
exciton creation operator are close to that of bosons 
and the corrections arise proportionally to the param- 
eter n(a B /£) 2 . 

Fig.Q]shows the behaviour of a n for different values of 
j e h interpolating from the bosonic case (jeh = —^\Jlelh) 
to the fermionic case (7^ = 0). The crossover from 
bosonic to fermionic limit can be clearly seen: for j e h 
close to — 2 A /7e7/ 1 , the curve behaves like y/n, the devia- 
tions from this exact bosonic result becoming more pro- 
nounced with increasing n. For 7^ close to —2^Jjejh, 
the curve initially behaves like y/n, the deviations from 
this exact bosonic result becoming more pronounced with 
increasing n. The curve is ultimately decreasing beyond 
a number of excitations which is smaller the greater the 
departure of j e h from —2^/7^7^. After the initial rise, as 
the overlap between electron and hole wavefunctions is 
small and bosonic behaviour is found, the decrease fol- 
lows as the density becomes so large that Pauli exclusion 
becomes significant. Then excitons cannot be considered 
as structure-less particles, and fermionic characteristics 
emerge. With 7^ going to 0, this behaviour is replaced 
by a monatonically decreasing a n , which means that it is 



n 




1 3 5 

n 



FIG. 4: (a) Matrix elements a n of the exciton creation op- 
erator X* calculated for n < 15 for various trial wave- 
functions. The top curve shows the limit of true bosons 
where ot n = \/n and the bottom curve the limit of true 
fermions where a n = S n ,i- Intermediate cases are obtained 
for values of y e h from —1.95^/7^7^ down to —0.2^/7^7^, inter- 
polating between the boson and fermion limit, (b) Magnified 
region close to the fermion limit. Values displayed are every- 
where given in units of ^y^jh. 



1 3 5 7 9 11 13 
n 

FIG. 5: Second order correlator 52(0) of the excitonic field 
for Fock states \n) given by Eqs. 115H16I as a function 
of 7 e h/^/7e7h which interpolates between the fermionic (—0.5, 
— 1 and —1.5) and bosonic (—1.9, —1.95, —1.99) limits. Up- 
per line corresponds to pure bosons (7^ = — 2^/7vyh). With 
increasing fermionic character the antibunching is very pro- 
nounced at high intensities. Close to the fermionic limit, Q2 
is always very small and vanishes quickly. 



dynamics so t — > 00. Eq. I|47|) reduces to 172(0) a quan- 
tity, which is one of great physical and experimental rel- 
evance. It will be sufficient for our description to con- 
sider Fock states of excitons only, although an exten- 
sion to other quantum states is straightforward. The 
matrix representation in the basis of states \n) reads 
X^X 2 = {a 2 _ 1 a 2 8i_j) Q <i_ J with a» = if i < 1, so that 
for a Fock state with n excitons, 



52(0) = (49) 



"harder and harder" to add excitons in the same state in 
the QD; the fermionic nature of excitons becomes more 
and more important. 

An important quantity for single mode particles, es- 
pecially in connection to their coherent features, is the 
(normalised) second order correlator g 2 which in our case 
reads 



92 (t, t) 



{X\t)X\t + T)X{t + T)X(t)) 

(A(t))(A(* + r)) 



(47) 



where p, is the exciton number operator which satisfies 



\i \n) 



(48) 



with |n) the bare exciton state with n electron-hole pairs, 
cf. (|15H16|) . At our energy diagram level, we can only 
compute zero-delay r = correlations and there is no 



Note that regardless of the model, 32(0) = for n = 1. 
General results are displayed on Fig.|3J They correspond 
to the exciton field which is the one of most interest to 
us, and could be probed in a two-photons correlation 
experiment with the light emitted directly by the exci- 
ton. We compare the result to the pure bosonic case 
where a n = \fn and therefore 52(0) = (n — l)/n so 
that 52 1 with increasing number of particles which 
expresses the similarity of an intense Fock state with a 
coherent state (especially regarding their fluctuations). 
In our case, however, the underlying fermionic structure 
results in an antibunching of excitons, i.e., the probabil- 
ity of finding two excitons at the same time is lowered at 
high exciton densities. Close to the fermionic limit, this 
antibunching is very pronounced and it is very unlikely 
to have more than one exciton in the system. 
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COUPLING TO A SINGLE RADIATION 
MODE OF A MICROCAVITY 



We now present the emission spectra of the coupled 
cavity-dot system when the exciton field is described by 
the creation operator X\ The procedure is straight- 
forward in principle and is a direct extension of the con- 
cepts discussed at length in the introduction. The hamil- 
tonian assumes the same form as previously, but now 
with X as wc defined it for the matter-field 6: 

H = hcj(a f a + fl) + hg(aX f + a' l X), (50) 

As (|50|l conserves the total number of excitations /j, = n+ 
m, it can be decoupled by decomposition of the identity 

as 



emission into a non-cavity mode. Then 



oo min(AT,/i) 

1 = Y^ \n, n - n) (n - n,\ 

p=0 n=0 



(51) 



where N is the smaller index n for which a n becomes 
zero. In the Gaussian wavefunction approximation with- 
out interaction, no a n ever becomes exactly zero, in 
which case N — > oo and the upper limit in the second 
sum of Eq. ifST]) is fi. 

Inserting l|51(l twice in (|50|l yields H = (B^L with 

Hn = hujfi + 



%KV (J> — n+ 15 n>m+ i + h.c.] 1 <„^„< min(A r iM ) (52) 
in the basis of bare photon-exciton states: 

W, t = {|0,M),|l,M-l),---, 

\imn(N, fi), fi - min(N, fj,)}} . (53) 

If a„ does not vanish, the basis has \i + 1 states in this 
manifold with last state \fJ-,0) having all photons trans- 
ferred in the excitonic field. In the case where a n van- 
ishes, the excitonic field saturates and the further exci- 
tations are constrained to remain in the photonic field. 

Following the nomenclature laid down in the introduc- 
tion, we write the dressed state \v, /j,} for the v tYl eigen- 
state of the manifold with /i = n + m excitations (n 
excitons + m photons) and = (n, v — n\v, ju) its de- 
composition on bare states of this manifold, i.e., 



min(N,fi) 



c£ M \n,fi 



(54) 



We compute the emission spectra corresponding to tran- 
sitions between multiplets, with matrix elements J en d = 
\{y',H — \\a\v,n)\ 2 for emission of a photon from the 
cavity and Ii at = \{v',fi — II^IKm)! 2 f° r direct exciton 



-^end — 



min(Af,;i-l) 
n=0 



'lat 



min(7V,/j) 

E (^ir 1 )*^ 



(55) 
(56) 



In cavity QED terminology, 7 cn d and ii a t correspond to 
end-emission and lateral-emission photo-detection re- 
spectively, while in luminescence of a microcavity, one 
observes the linear combination of the two contributions 
simultaneously. I\ a t reflects most the behaviour of the ex- 
citonic field. To separate the two, one can make use of the 
scattering geometry; in a pillar structure, for instance, 
the cavity photon emission is predominantly through the 
end mirrors, so I\ at could be detected on the edge of the 
structure. These measurements display the most inter- 
esting features and we focus on them. 

Fig. HJ1 shows the calculated emission spectra 7i at of a 
QD embedded in a cavity. All the spectra are broad- 
ened by convolution with a Lorentzian of width 7 = 
0.2<7. Figs, (a), (b) and (c) are, respectively, the re- 
sults close to the bosonic limit (7^ = — 1.99^/7 e 7^), 
in between {^ e h = —^Jlelh) and close to the fermionic 
limit (7c/! = — 0.05^/7e7^). Each curve is labelled with 
the number of excitations in the manifold, reflecting the 
intensity of the pumping field. Fig. EI a) shows a pro- 
nounced Rabi doublet in the case n = 1 (vacuum field 
Rabi splitting) in accordance with the general theory set 
out in the introduction. Higher manifolds reveal non- 
bosonic behaviour, with a reduced Rabi splitting doublet 
when n = 2, and the onset of a multiplct structure when 
n = 5, where small peaks appear at E w h(ui ± 2g). The 
increase of dot confinement in the case j e h = — ^/"felh 
in Fig. Ei|b), makes deviations from bosonic behaviour 
more pronounced, so even for n = 2 a multiplet struc- 
ture is observable. For higher manifold numbers, more 
lines appear, and a decrease of the splitting of the central 
Rabi doublet occurs. Further strengthening the confine- 
ment to 7e/j = — 0.05^/7 e 7? l makes the manifestations 
of fermionic behaviour very clear; for n — 2, a quadru- 
plet structure is seen, with the central two peaks almost 
merged. The situation is similar to the Mollow triplet of 
the exact fermionic limit, even more so for an increase of 
manifold number from n = 2 to n = 10, where the sepa- 
ration of the side peaks grows and the central ones merge 
further. It can be seen from Fig. El that the decrease of 
the QD size and the corresponding change of the exciton 
quantisation regime (with 7^ ranging from —2 y /j e ^f] l to 
0) is manifest in the emission spectra as a transition from 
a Rabi doublet to a Mollow triplet. The increase of the 
excitation power (at fixed QD size) leads qualitatively to 
the same effect; fermionic behaviour becomes more pro- 
nounced with the increase of the excitation power and 
decrease of the QD size. 

Having demonstrated how our model interpolates be- 
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FIG. 6: (Color online) Spectra for various intensities of the 
light field at specific values of 7^: (a) close to the bosonic 
limit with j e h = — 1.9-v/7e7h for n = 1 (solid red), 2 (dashed 
blue) and 7 (dotted black) featuring a broadened and red- 
shifted Rabi doublet as the intensity increases and the onset 
of a multiplet structure, (b) intermediate case: -y e h = — ^/7e7h 
for n — 2 (solid black) and 10 (dashed blue) demonstrating a 
complicated multiplet structure and (c) close to the fermionic 
case with 7^ = — 0.05^/ r y e yh for n = 2 (solid black) and 10 
(dashed blue) featuring quadruplet structure going towards 
Mollow triplet at high intensities. 



tween Fermi and Bose statistics, we now investigate in 
further detail the intermediate regime, by following the 
evolution of a single manifold between the two limits. 
In Fig. the eigenvalues are displayed for the 7th man- 
ifold, also shown on Fig. ^ f° r bosons and fermions. 
Eigenvalues are plotted as a function of y e h running 
from —2^/7^7^, which recovers the left hand side of Fig. ^ 
to 0, which recovers the right hand side. Other mani- 
folds behave in qualitatively the same way. In all cases, 
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FIG. 7: Energy diagram of the dressed excitons over the in- 
terval 7 e fe £ [— 2^/7 e 7h, 0] for the 7th manifold. The extremal 
values recover the diagrams of Fig.Qwith n + 1 equally spaced 
energy levels in nth dressed boson manifold (far left) and 
twofold energy diagram with square root splitting for dressed 
fermions (far right). The two outer boson energies connect 
smoothly to the two fermion energies while the n — 1 other 
one degenerate into a central line which disappears right when 
the system hits the Fermi limit (this point, 7^ = 0, is shown 
by circle). 



the n + 1 equally spaced energies of the dressed bosons 
link to the two energies of the dressed fermions as follows: 
the upper energy /ihuj + \1T1g of the Bose limit links to 
the upper energy fihtu + y/JIhg of the Fermi limit, and the 
symmetric behaviour occurs with the lower limit, link- 
ing fihuj — fihg with [ifiLo — yJJLTig. More interestingly, the 
intermediate energies degenerate from the equal spacing 
of the Bose limit towards [ihw near the Fermi limit, with 
a discontinuity when 0L2 becomes zero and only two eigen- 
cnergies remain. Physically, this discontinuity arises be- 
cause for any 7,,^ 7^ 0, formally, an infinite number of 
excitations can be fitted in the QD, i.e. a n never be- 
comes exactly zero, see Fig. 0] With j e h approaching 
zero these modes couple to light more and more weakly, 
thus the energies of the dressed states become close to the 
bare energy (ifitu. When 7^ = the fermionic regime is 
recovered and only two dressed eigenstates with energies 
fj,hui ± ^/JIh,g survive. 

Although the eigenvalue structure of the dressed ex- 
citons in the general case is simple and uniform, as is 
seen in Fig. [7] the multiplet structures which arises from 
it is, as has been seen in the spectra previously in- 
vestigated, rich and varied. In the general case where 
— 2^/7 e 7/i < 7 e h < 0, two adjacent manifolds contain re- 
spectively n + 1 and n levels, transitions between any pair 
of which are possible. One may expect n(n + 1) lines in 
the emission spectra with a range of intensities and posi- 
tions. The phenomenological broadening we have intro- 
duced leads to the decrease of the number of resolvable 
lines. In all cases the n = 1 to n = transition pro- 
vides only two lines, with constant energies throughout, 



12 




Boson Fermion 



FIG. 8: Transition from Fermi to Bose limits as observed in 
the optical emission spectra, (a) Superpositions of spectra 
for "feh = — 2 x /7 e 7h in front of the figure to 7 e h = at the 
back, recovering respectively the Rabi doublet and Mollow 
triplet. In the intermediate region, intricate and rich patterns 
of peaks appear, split, merge, or disappear, (b) Same data as 
in (a), but as a density plot on logarithmic scale, to discrim- 
inate the peaks, their positions and splitting as well as their 
relative intensities. All spectra are normalized to maximum 
intensity equal to 1. 



As a final comment, we note that if Eq. i|46[) was to 
hold for all n, it would yield, apart from a renormaliza- 
tion by y/N, the Dicke models, which has been widely 
used to describe various strong light-matter coupling phe- 
nomena-2iii. In our model it is recovered when Pm_= 
l/iV( m-1 ), a case that we have investigated in Ref. [421 ]. 

At the heart of Dicke model lies the creation oper- 
ator J + for an excitation of the "matter field" which 
distributes the excitation throughout the assembly of A 
identical two-levels systems described by fermion opera- 
tors <7i, so that in Eq. @ maps to J + with 

N 

J+ = 5>? (57) 

i=l 

One checks readily that J + and J_ = j| thus de- 
fined obey an angular momentum algebra with magni- 
tude N(N + 1) (and maximum z projection of J z equal 
to A). In this case the Rabi doublet arises in the 
limit where the total number of excitations [i (shared 
between the light and the matter field) is much less 
than the number of atoms, /i <C A, in which case the 
usual commutation relation [J_,J + ] = — 2J Z becomes 
\J-/\f~N, Jj r /\/~N] ~ 1, which is the commutation for a 
bosonic field. This comes from the expression of a Dicke 
state with \i excitations shared by A~ atoms given as the 
angular momentum state \—N/2 + /i). Therefore the an- 
nihilation/creation operators J_, J+ for one excitation 
shared by N atoms appear in this limit like renormalised 
bose operators \ Na, \/~Na\ resulting in a Rabi doublet 
of splitting 2hg\/N. Such a situation corresponds, e.g., 
to an array of small QDs inside a microcavity such that in 
each dot electron and hole are quantized separately, while 
our model describes a single QD which can accommodate 
several excitons. The corresponding emission spectra are 
close to those obtained here below the saturation limit 
/j <C N, while the nonlinear regime JV > 1, /j » 1 
has peculiar behaviour, featuring non-lorentzian emission 
line shapes and a non-trivial multiplct structure, like the 
"Dicke fork" M 



VI. CONCLUSIONS AND PROSPECTS 



because with only one exciton present the question of 
bosonic or fermionic behaviour is irrelevant. To access 
the consequences of exciton statistics, one must, unsur- 
prisingly, reach higher manifolds. 

In Fig. [S] we show the evolution of the spectra between 
the boson and fermion limits. They two panels are dif- 
ferent projections of the same data, namely the multiplct 
structure as a function of j e h € [—2^/7^7^", 0] for n = 12. 
The Rabi doublet is seen to evolve into a Mollow triplet 
going through a complex and intertwined set of peaks 
whose splitting and relative heights vary with the value 
of 7 e ^ considered. The spectrum obtained is therefore a 
direct probe of the underlying exciton quantum statistics. 



The spectrum of light emitted by quantum dot exci- 
tons in leaky modes of a microcavity — which could be 
typically lateral emission — is a signature of the quantum 
statistics of excitons. A multiplct structure has been the- 
oretically predicted with various features which can help 
identify the exciton field statistics. We provided the for- 
malism to obtain the spectra expected for a general dot 
in various geometries based on the form of the single exci- 
ton wavefunction. We investigated a generic case analyt- 
ically through a Gaussian approximation. The richness 
and specificity of the resulting spectra provides a means 
to determine, through the peak splittings and strength 
ratios, the parameter 7^ that measures how close exci- 
tons are to ideal fermions or bosons. Although heavy nu- 
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merical computations arc required for realistic structures, 
physical sense motivates that in small dots a Fermi-like 
behaviour of excitations with separately quantised elec- 
trons and holes is expected, while in larger dots behaviour 
should converge towards the Bose limit. These trends 
should be observable in the nonlinear regime (where more 
than one exciton interacts at a single time with the ra- 
diation mode): depending on whether Rabi splitting is 
found, or if a Mollow triplet or a more complicated mul- 
tiplet structure arises, one will be able to characterize the 
underlying structure of the exciton held. 

We have left the discussion of the interaction between 
excitons out of the scope of the present paper consid- 
ering the effects which arise solely from Pauli exclusion 
principle. The effect of interactions, through screening of 
Coulomb potential by the electron- hole pairs, has been 
discussed in Refs. [37]E3- These papers demonstrate 
that the interactions will make fermionic behaviour more 
pronounced, since presence of other excitons screens the 
Coulomb interaction which binds electron and hole and 
leads to the increase of the Bohr radius. The question 
whether interactions will predominate over Pauli exclu- 
sion lies beyond the scope of this paper and we post- 



pone it for a future work. Our preliminary estimations 
show that even at moderate excitonic density the effect of 
screening is small and does not lead to strong qualitative 
deviations^ 

To summarise, we have studied the effect of Pauli 
exclusion on the optical emission spectra of microcavi- 
ties with embedded QDs in the strong coupling regime. 
We derived general expressions for the exciton creation 
operator which allow systematic computation of the 
light-matter coupling. The crossover between bosonic 
behaviour — observed in large QDs as Rabi doublet — to 
the fermionic behaviour — observed in small QDs as Mol- 
low triplet — has been demonstrated. 
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